Dzyaloshinsky-Moriya antisymmetric exchange coupling in cuprates: Oxygen effects 



o 
o 

(N 

> 

O 

ON 



CO 



I 

C 

O 
O 



> ■ 

■ 

cn : 
in . 
o ■ 

: 

o . 



c : 
o 

o ■ 



X 



A.S. Moskvin 1 

J f/raZ SZate University, 620083 Ekaterinburg, Russia 
(Dated: February 4, 2008) 

We revisit a problem of Dzyaloshinsky-Moriya antisymmetric exchange coupling for a single bond 
in cuprates specifying the local spin-orbital contributions to Dzyaloshinsky vector focusing on the 
oxygen term. The Dzyaloshinsky vector and respective weak ferromagnetic moment is shown to 
be a superposition of comparable and, sometimes, competing local Cu and O contributions. The 
intermediate oxygen 17 O Knight shift is shown to be an effective tool to inspect the effects of 
Dzyaloshinsky-Moriya coupling in an external magnetic field. We predict the effect of strong oxygen 
weak antiferromagnetism in edge-shared CuCh chains due to uncompensated oxygen Dzyaloshinsky 
vectors. Finally, we revisit the effects of symmetric spin anisotropy, in particular, those directly 
induced by Dzyaloshinsky-Moriya coupling. 



I. INTRODUCTION 



Fifthy years ago Borovik-Romanov and Orlova 1 
have discovered the phenomenon of weak ferromag- 
netism which origin was shortly after— bound up with 
exchange-relativistic effect mainly with antisymmetric 
exchange coupling. Starting from pioneer papers by 
Dzyaloshinsky 2 and Moriya^ the Dzyaloshinsky-Moriya 
(DM) antisymmetric exchange coupling was extensively 
investigated in 60-80ths in connection with weak ferro- 
magnetism focusing on hematite cv-Fe203 and orthofcr- 
rites RFe03>i The stimulus to a renewed interest to the 
subject was given by the cuprate problem, in particu- 
lar, by the weak ferromagnetism observed in La2CuO.$ 5 
and many other interesting effects for the DM systems, 
in particular, the "field-induced gap" phenomena— At 
variance with typical 3D systems such as orthoferrites, 
cuprates are characterised by a low-dimensionality, large 
diversity of Cu-O-Cu bonds including corner- and edge- 
sharing, different ladder configurations, strong quantum 
effects for s — 1/2 Cu 2+ centers, and a particularly strong 
Cu-0 covalency resulting in a comparable magnitude of 
hole charge/spin densities on copper and oxygen sites. 
Several groups (see, e.g., RefsQIHa) developed the micro- 
scopic model approach by Moriya for different ID and 2D 
cuprates, making use of different perturbation schemes, 
different types of low-symmetry crystalline field, differ- 
ent approaches to intra-atomic electron-electron repul- 
sion. However, despite a rather large number of publica- 
tions and hot debates (see, e.g., Refllfl) the problem of 
exchange-relativistic effects, that is of antisymmetric ex- 
change and related problem of spin anisotropy in cuprates 
remains to be open (see, e.g., Refs lllll^ for recent exper- 
imental data and discussion). Common shortcomings of 
current approaches to DM coupling in 3d oxides concern 
a problem of allocation of Dzyaloshinsky vector and re- 
spective "weak" (anti)ferromagnetic moments, and full 
neglect of spin-orbital effects for "nonmagnetic" oxygen 
2 ~ ions, which are usually believed to play only indi- 
rect intervening role. From the other hand, the oxygen 
17 NMR-NQR studies of weak ferromagnet L^CuO^ 
seem to evidence unconventional local oxygen "weak- 



ferromagnetic" polarization which origin cannot be ex- 
plained in frames of current models. All this stimulated 
the critical revisit of many old approaches to the spin- 
orbital effects in 3d oxides, starting from the choice of 
proper perturbation scheme and effective spin Hamilto- 
nian model, implied usually only indirect intervening role 
played by "nonmagnetic" oxygen 2 ~ ions. 

In this paper we revisit a problem of DM antisymmet- 
ric exchange coupling for a single bond in cuprates spec- 
ifying the local spin-orbital contributions to Dzyaloshin- 
sky vector focusing on the oxygen term. In Sec. II we 
present a short overview of the spin-Hamiltonian of a 
typical three-center (Cu-O-Cu) two-hole system. Micro- 
scopic theory of DM coupling is presented in Sec. III. The 
Dzyaloshinsky vector is shown to be a superposition of 
comparable and, sometimes, competing local Cu and O 
contributions. In Sec. IV we examine a response of DM 
coupled Cui-0-Cu2 bond on an uniform and staggered 
external fields, and demonstrate some unusual manifes- 
tations of local oxygen contribution to DM coupling. The 
intermediate oxygen 17 O Knight shift is shown to be an 
effective tool to inspect the effects of DM coupling in an 
external magnetic field. In Sec.V we revisit a problem of 
symmetric spin anisotropy with inclusion of local oxygen 
spin-orbital contributions. 




FIG. 1: Geometry of the three-center (Cu-O-Cu) two-hole 
system with ground Cu 2>A x 2_ y 2 states. 
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II. SPIN-HAMILTONIAN 

For illustration, below we address a typical for cuprates 
the three-center (Cui-0-Cu 2 ) two-hole system with 
tetragonal Cu on-site symmetry and ground Cu 3d x 2_ y 2 
states (see Fig. 1) which conventional bilinear spin 
Hamiltonian is written in terms of copper spins as fol- 
lows 

H s (12) = J u (si ■ s 2 ) + D 12 • [sj x s 2 ] + sxKu §2 , (1) 

where J12 > is an exchange integral, D12 is the 

Dzyaloshinsky vector, K12 is a symmetric second-rank 
tensor of the anisotropy constants. In contrast with 

Ji2,Ki2, the Dzyaloshinsky vector D 12 is antisymmet- 
ric with regard to the site permutation: D12 = — D21. 

I 



which implies not only both copper and oxygen hole lo- 
cation, but allows to account for purely oxygen two-hole 
configurations. Moreover, such a form allows us to neatly 
separate both one-center and two-center effects. Two- 
hole spin Hamiltonian can be projected onto three-center 
states incorporating the Cu-0 charge transfer effects. 

For a composite two s = 1/2 spins system one should 
consider three types of the vector order parameters: 

S = si+s 2 ; V = si-s 2 ;f = 2[sixs 2 ] (3) 

with a kinematic constraint: 

S 2 +V 2 = 31; (S-V) = 0; (f -V) = 6i; [TxV] = S. (4) 

In a sense the V operator describes the effect of local 
antiferromagnetic order, while T operator may be asso- 
ciated with a vector chirality^ In recent years, phases 
with broken vector chirality in frustrated quantum spin 
chains have attracted considerable interest. Such phases 
are characterized by nonzero long-range correlations of 
the vector order parameter (T) . Interestingly that a chi- 
rally ordered phase can manifest itself as a "nonmag- 
netic" one, with (S) = (V) = 0. 

Both T and V operators change the spin multiplicity 
with matrix elements 

(00|f m |ln) = -<ln|f m |00) = iS mn ; 

(00|V;„|ln) = (ln|t/ m |00) = S mn , (5) 

where we made use of Cartesian basis for S = 1. The 
eigenstates of the operators V n and T n with nonzero 



Hereafter we will denote J i2 = J, K12 = K, Di 2 = D, re- 
spectively It should be noted that making use of effective 
spin Hamiltonian Q implies a removal of orbital degree 
of freedom that calls for a caution with DM coupling as 
it changes both a spin multiplicity, and an orbital state. 

It is clear that the applicability of such an operator as 
-ff s (12) to describe all the "oxygen" effects is extremely 
limited. Moreover, the question arises in what concerns 
the composite structure and spatial distribution of what 
that be termed as the Dzyaloshinsky vector density. Usu- 
ally this vector is assumed to be located on the bond 
connecting spins 1 and 2. 

Strictly speaking, spin Hamiltonian H s (12) can be 
viewed as a result of the projection onto the purely ionic 
Cu 2+ (3d x 2_ y 2)-0 2 ~ (2p 6 )-Cu2 + (M x 2_ y 2) ground state of 
the two-hole spin Hamiltonian 



(2) 

I 

eigenvalues ±1 form Neel doublets "^(|00) ± \ln)) and 
DM doublets ^(|00) ± i\ln)), respectively. The Neel 
doublets correspond to classical collinear antiferromag- 
netic spin configurations, while the DM doublets corre- 
spond to quantum spin configurations which sometimes 
are associated with a rectangular 90° spin ordering in the 
plane orthogonal to the Dzyaloshinsky vector. 

It should be noted that the both above spin Hamil- 
tonians can be reduced to within a constant to a spin 
operator acting in a total spin space 

H s = ij(S 2 -V 2 ) + i(D-t) + isK S S-ivK y V. (6) 

For simple dipole-like two-ion anisotropy as in Exp. (1) 

K = K = K, though in general these tensorial pa- 
rameters can differ from each other. Making use of the 
anticommutator relations 

S 3 } + {% Vj} = 2%; {Vi, Vj} = {futj}, (7) 

we conclude that the effective operator of symmetric 
anisotropy can be equivalently expressed in terms of sym- 
metric products {§i,Sj}, {Vi,Vj}, or {Tj,Tj}. 

The most general form of the spin Hamiltonian © 
does not discriminate against copper or oxygen contri- 
bution and can be used to properly account for oxygen 

effects. As we will see below the D and K parameters 
allow the correct separation of local copper and oxygen 
contributions. 

It should be noted that the total spin representation 
and overall quantum approach is especially efficient to de- 
scribe antisymmetric exchange in Cu-Cu dimer systems 



i<j 



i<j 
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(see, e.g. ReflTa and references therein). Classical ap- 
proach to s = 1/2 spin-Hamiltonian should be applied 
with a caution, particularly for ID systems. 

Before going to microscopic analysis we should note 
that the interaction of our three-center system with ex- 
ternal spins and/or fields H ext is usually addressed by 
introducing only two types of effective external fields: 
the conventional Zeeman-like uniform field and uncon- 
ventional Neel-like staggered field, so that H ex t reads as 
follows 



H ext = -(h s ■ S) - (h v -V). 



(8) 



It should be noted that an ideal Neel state is attainable 
only in the limit of infinitely large staggered field, there- 
fore for a finite staggered field h v | n the ground state 
is a superposition of a spin singlet and a Neel state, 



if? = cosa|00) + sinajln), tan 2a 



2hy 
j 



which composition reflects the role of quantum effects. 
For instance, in a Heisenberg spin 1/2 chain with nn ex- 
change the maximal value of staggered field h v = J/2 
hence the \& function strongly differs from that of Neel 
state ((V n ) — sin 2a = -^=), and quantum mechanical av- 
erage for a single spin (s z ) < ^ sin7r/4 = ■ | w 0.71 • | 

deviates strongly from classical value |. It should be 
noted that for the isolated antiferromagnetically coupled 
spin pair the zero-temperature uniform spin susceptibil- 
ity turns into zero: x S = 0? while for the staggered spin 
susceptibility we obtain \ V — 2/ J. 



III. MICROSCOPIC THEORY OF DM 
COUPLING IN CUPRATES 

A. Preliminaries 

For the microscopic expression for Dzyaloshinsky vec- 
tor to derive Moriya— made use of the Anderson's formal- 
ism of superexchange interaction^ with two main con- 
tributions of so called kinetic and potential exchange, 
respectively. Then he took into account the spin-orbital 
corrections to the effective d-d transfer integral and po- 
tential exchange. Such an approach seems to be im- 
proper to account fo r pur ely oxygen effects. In later pa- 
pers (see, e.g. Refs l8ll7T) the authors made use of the 
Moriya scheme to account for spin-orbital corrections to 
p-d transfer integral, however, without any analysis of 
oxygen contribution. It is worth noting that in both in- 
stances the spin-orbital renormalization of a single hole 
transfer integral leads immediately to a lot of problems 
with correct responsiveness of the on-site Coulomb hole- 
hole correlation effects. Anyway the effective DM spin- 
Hamiltonian evolves from the high-order perturbation ef- 
fects that makes its analysis rather involved and leads to 
many misleading conclusions. 



At variance with the Moriya approach we start with 
the construction of spin-singlet and spin-triplet wave 
functions for our three-center two-hole system taking ac- 
count of the p-d hopping, on-site hole- hole repulsion, and 
crystal field effects for excited configurations {n} (Oil, 
110, 020, 200, 002) with different hole occupation of Cui, 
O, and Cu2 sites, respectively. The p-d hopping for Cu-0 
bond implies a conventional Hamiltonian 



H p d — 



(9) 



a/3 



where f>\ creates a hole in the a state on the oxygen 
site, while dp annihilates a hole in the (3 state on the 
copper site; t pa df3 is a pd-transfer integral {tp x d x2 _ i = 

For basic 101 configuration with two d x 2_ y 2 holes lo- 
calized on its parent sites we arrive at a perturbed wave 
function as follows 



101;SM 



= <f> 



Wl-SM 



{n}T 



C{n}( 



2S+1 



r)$ 



{n};TSM 



> (io) 



where the summation runs both on different config- 
urations and different orbital T states. It is worth 
noting that the probability amplitudes C{oii}, C{no} oc 
£prf,c { 20o}7C{020}:C{oo2} oc t 2 pd . To account for orbital ef- 
fects for Cui.2 3d holes and the covalency induced mixing 
of different orbital states for 101 configuration we should 
introduce an effective exchange Hamiltonian 



1 



J(a^)dl QjU 4^'^2 7M ^i<5M'+' l - c - ( n ) 



Here d\ ali creates a hole in the ath 3d orbital on Cui site 
with spin projection [i. Exchange Hamiltonian Ijlll) in- 
volves both spinless and spin-dependent terms, however, 
it preserves the spin multiplicity of Cui-0-Cu2 system. 
Exchange parameters J(af3jS) are of the order of t pd . 

Then we introduce a standard effective spin Hamil- 
tonian operating in a fourfold spin-degenerated space 
of basic 101 configuration with two d x 2_ y 2 holes and 
can straightforwardly calculate the singlet-triplet sep- 
aration to find effective exchange integral Jyi — 
J(d x 2_ y 2d x 2_ y2 d x 2_ y 2d x 2_ y 2) 1 calculate the singlet- 
triplet mixing due to three local spin-orbital terms, 
V so (Cui), V so (0), V so (Cu 2 ), respectively, to find the 
local contributions to Dzyaloshinsky vector: 



d = d« + d(°)+d( 2 ). 

Local spin-orbital coupling is taken as follows: 



(12) 



V so = J2^i(h ■ Si) = ^[(li + la) • S + (h - h) ■ V] 



V 



(13) 
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with a single particle constant > for electrons 
and < for holes. We make use of orbital 

matrix elements: for Cu 3d holes (d x a_ y a \l x \d yz ) = 
(d x 2_ y 2\l y \d xz ) = i, (d x 2_ v 2\l z \d xy ) = -2i, (i\lj\k) = 
-ie ljk with Cu 3d^ = |l), S& xz =\2) 3d a:v =|3), and for 
2p holes (pi\lj\pk) — itijk- Free cuprous Cu 2+ ion is de- 
scribed by a large spin-orbital coupling with \^zd\ = 0.1 
eV (see, e.g., Refll8r). though its value may be signifi- 
cantly reduced in oxides. Information regarding the ^2p 
value for the oxygen O 2- ion in oxides is scant if any. 
Usually one considers the spin-orbital coupling on the 
oxygen to be much smaller than that on the copper, 
and therefore may be neglect edi^SiSl However, even for 
a free oxygen atom the electron spin orbital coupling 
turns out to reach of appreciable magnitude: ^2p — 0.02 
eVj^i while for the oxygen 2 ~ ion in oxides one ex- 
pects the visible enhancement of spin-orbital coupling 
due to a larger compactness of 2p wave function^ If to 
account for cx (r~ 3 ) n ; and compare these quantities 
for the copper ((r -3 ^ « 6 — 8 a.u.— ) and the oxygen 
{( r ~ /2p ~ 4 a.u^ 3 * 2 ^) we arrive at a maximum factor 
two difference in £3^ and ^2p- 

Hereafter we assume a tetragonal symmetry at Cu sites 
with local coordinate systems as shown in Fig.l. The 
global xyz coordinate system is chosen so as CU1-O-CU2 
plane coincides with xy-plane, x-axis is directed along 



In all the instances the nonzero contribution to the local 
Dzyaloshinsky vector is determined solely by the spin- 
orbital singlet-triplet mixing for the on-site 200, 020, 
002 and two-site 110, 011 two-hole configurations, re- 
spectively. For on-site two-hole configurations we have 

D (200) = D (l) ) D (020) = D (0) j D (002) = D (2)_ 
I 



1 ' should be noted that at variance with the origi- 
nal Moriya approach^ both spinless and spin-dependent 
parts of exchange Hamiltonian contribute additively and 
comparably to DM coupling, if one takes account of the 
same magnitude and opposite sign of the singlet-singlet 



CU1-CU2 bond (see Fig.l). In such a case the basic unit 
vectors x, y, z can be written in local systems of Cui and 
CU2 sites as follows: 

x = (sin - , — cos - cos 01 , — cos - sin di); 



y = (cos — , sin — cos 8%, sin — sin<5i); z = (0, sin 8x, cos 8±) 

for Cui, while for Cu2 site 9,8\ should be replaced by 
—6,82, respectively. 

The exchange integral can be written as follows: 

J = E[l c w( 3r )| 2i?3 r({n})-|c w ( 1 r)| 2 £; lr ({n})]. 
{n},r 

(14) 

As concerns the DM interaction we deal with two com- 
peting contributions. The first is derived as a first order 
contribution which does not take account of Cui^ 3d- 
orbital fluctuations for a ground state 101 configuration. 
Projecting the spin-orbital coupling (|13|l onto states I|1U|) 
we see that A v ■ V term is equivalent to purely spin DM 
coupling with local contributions to Dzyaloshinsky vector 



(15) 

I 

The second contribution, associated with Cui^ 3d- 
orbital fluctuations within a ground state 101 configura- 
tion, is more familiar one and evolves from a second order 
combined effect of Cui^ spin-orbital V so (Cu\,2) and ef- 
fective orbitally anisotropic CU1-CU2 exchange coupling 



(16) 



and triplet-triplet exchange matrix elements on the one 
hand and orbital antisymmetry of spin-orbital matrix el- 
ements on the other hand. 

It is easy to see that the contributions of 002 and 200 
configurations to Dzyaloshinsky vector bear a similarity 



J 



Dt } = -2i(00\V so (m)\li) = -2i £ ^ n} ( 1 r 1 ) C{ „ } ( 3 r 2 )($ {ll}; r 1 oo|Ar |* {n}; r a ii> 

{n}ri,r a 



9 /nnlT/ < \\m\ «V ^ 101 > ; r B 00|AK|{lQl}; riz)({101}; riz|g e:c |{101}; r f lz) 
MMWso{m)\U) = -2.^ ift rt ({l01})-iM{101}) 



- 2 *E 

r 



({101}; r s 00|g ex |{101}; TOO) ({101}; r00| Af |{101}; T t li) 
^r,({101})-^ r ({101}) 

I 
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to the respective second type (oc V so H ex ) contributions, 
however, in the former we deal with spin-orbital coupling 
for two-hole Cui^ configurations, while in the latter with 
that of one-hole Cui ; 2 configurations. 



B. Copper contribution 

First we address a relatively simple instance of strong 
rhombic crystal field for intermediate oxygen ion with the 
crystal field axes oriented along global coordinate x, y, z- 
axes, respectively. It is worth noting that in such a case 
the oxygen O 2p z orbital does not play an active role 
both in symmetric and antisymmetric (DM) exchange 



J 



interaction as well as Cu 3d yz orbital appears to be in- 
active in DM interaction due to a zero overlap/transfer 
with oxygen O 2p orbitals. 

For illustration, hereafter we address the first contribu- 
tion H15|) of two-hole on-site 200, 002 d x2 _ y2 , d x 2_ y 2d xy , 
and d x i_yid xz configurations, which do covalently mix 
with ground state configuration. Calculating the singlet- 
triplet mixing matrix elements in the global coordi- 
nate system we find all the components of the local 
Dzyaloshinsky vectors. The Cui contribution turns out 
to be nonzero only for 200 configuration, and may be 
written as a sum of several terms. First we present the 
contribution of singlet {d x i _ y i) x A\ g state: 



^ X 



Q 

-2i(00|V r 8O (£7ui)|lx) = \ / 26dC2oo( 1 A lg )[c 20 o( 3 E g ) cos ft - 2c 2ao ( 3 A 2g ) sin ft] cos -; 



V 



-2i{Q0\V so (C Ul )\ly) = -V2&dC 2 oo{ 1 A lg )[c 20 o( 3 'E g ) cos ft - 2c 20 o( 3 A 2g ] 



sin ft] sin — ; 



D« = -2i(00|V r . o (Cu 1 )|l«) = -V2& d C2oo{ 1 A lg )[c 2 oo( 3 E g ) sin 6 X - 2c 200 ( 3 A 2g ) cos ft], 

I 



(17) 



where 

C20o( 1 -4i 9 ) 



2V2 pda E lAlg 



V3. 



sm 



cos 



2 e 



C20o( ' A 2g ) — —tpda-tpdTT 



£1.3 Aa 



C20o( 1,3 ^g) 



V3 1 

4 Ei,s E 



1 1 1 • 

1 ) sin cos ft, 

e x e. 



sin 9 sin ft, 



are probability amplitudes for singlet (d^ 2 -j^) 1 Ai g and 
singlet/triplet {d x 2^ y 2d xy ) 1 > 3 A 2g , (d x 2_ y 2d X z) 1 ' 3 E g 200 
configurations in the ground state wave function, re- 
spectively. Ei a = A + 4£> + 3C is the energy of 



the two-hole copper singlet with d 2 _ 2 configuration, 
£m 2o = e xy + A + 4B + 2C, Es A * g = e xy + A + AB, 



E 



A + B + 2C, E 3E 



A - 5B 



are the energies of the two-hole copper terms with 
d x 2_ y 2d xy and d x 2_ y 2d xz configurations, respectively, A, 
B, and C are Racah parameters. Taking into account 

that C 002 ( 1 A lg ) = C 20 o( 1 ^lff), C002( 3 ^2g) = C 2 Oo( 3 ^2g), 

coo2( 3 E 2g ) — c 2 oo( 3 E g )m& we see that the CU2 contribu- 
tion to the Dzyaloshinsky vector can be obtained from 

Exps. ifTTj) . if 9, ft replace by —0,S 2 , respectively. In- 

(1 2) 

terestingly to note that D x ^ y oc sin9 sin2ft,2- Both 
collinear (9 = n) and rectangular (9 = n/2) superex- 
change geometries appear to be unfavorable for copper 



contribution to antisymmetric exchange, though in the 
latter the result depends strongly on the relation between 
the energies of oxygen O 2p2; and O 2pj, orbitals. 

Contribution of singlet {d x 2_ y 2d xy ) 1 A 2g and 
{d x 2_ y 2d xz ) l E g states to the Dzyaloshinsky vector 
yields 

<£>=<&> S m * 4 1 )=dW C os^, d«=0, 

where rf (1) = ^dic^ 1 A 2g )c 2m { 3 E g ) - 

C20o( 1 ^£/)c20o( 3 ^2 g ))- Here we deal with a vector directed 
along Cui-0 bond which modulus dP^ oc sin 2 9 sin 2ft is 
determined by a partial cancellation of two terms. 

It is easy to see that the copper V so (l) contribution to 
the Dzyaloshinsky vector for two-site 110 and 011 con- 
figurations is determined by a dp-exchange. 



C. Oxygen contribution 

In frames of the same assumption regarding the orien- 
tation of rhombic crystal field axes for intermediate oxy- 
gen ion the local oxygen contribution to the Dzyaloshin- 
sky vector for one-site 020 configuration appears to be 
oriented along local O z axis and may be written as fol- 
lows 



£>£°> = -2i(00\V so (O)\lz) = V2t; 2p c t (p xPy )[c( P l) 



(18) 
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where c(p 2 ) = co2o(p 2 ) with 
/ 2\ _ 3 2 sm 2 /• 2\ _ 3 2 cos2 2 

. 3 2 /I 1 \ sin6» 

= g W + Tj em (19) 

are probability amplitudes for singlet p 2 ,p 2 and triplet 
PxPy 020 configurations in the ground state wave func- 

I 



ft is worth noting that D^ ' does not depend on the Si, S 2 
angles. The Dq{9) dependence is expected to be rather 
smooth without any singularities for collinear and rect- 
angular superexchange geometries. 



where c s ,t(dp) = cno(dp),c Stt (pd) = cou(dp) are proba- 
bility amplitudes for different singlet (c s ) and triplet (q) 
ffO (d x 2_ y 2p xy ) and Oil (p x ,yd x 2^ y 2) configurations in 
the ground state wave function, respectively. The ener- 
gies E s ^{dp X: y) are those for singlet and triplet states of 
dp x .y configurations, respectively: E Si t(dp x , y ) = tx,y + 



where we take account only of dpa exchange (Jdpa oc 

Thus, the total Dzyaloshinsky vector D is a superpo- 
sition of three contributions D = + D<W + at- 
tached to the respective sites. In general, all the vectors 
can be oriented differently. Their magnitudes strongly 



tion, respectively; E s {p xy ) = 2e x<y + F + ^F 2 , 
E t (p x p y ) = e x + c y + Fq — ^F 2 are the energies of the 
oxygen two-hole singlet (s) and triplet (t) configurations 
p x ,p y and p x p y , respectively, Fq and F 2 are Slater inte- 
grals. This vector can be written as follows^ 

D<°> =D {6)[t 1 xr 4 ], (20) 

where ri j2 are unit radius-vectors along Cui-O, Cu2-0 
bonds, respectively, and 



cos 2 | sin 2 | 

^M) ~ ' (21) 

I 

The local oxygen contribution to the Dzyaloshinsky 
vector for two-site 110 and Oil configurations appears 
to be oriented along local O z axis as well and may be 
written as follows 



(22) 



(23) 

I 

Kd px , v ± Id P x,y, where Kd px , v and Id P x, v are Coulomb 
and exchange dp-integrals, respectively. It is easy to see 
that the nonzero contribution to Dzyaloshinsky vector is 
determined by a direct dp-exchange and may be written 
similarly to l|2U|l with 



(24) 

I 

depend on the Cui-0-Cu2 bond geometry and crystal 
field effects, however, comparative analysis of Exps. i|17|) 
and (|21ll given estimations for different parameters typi- 
cal for cupratesSi {t p da ~ 1-5 eV, t v d^ ~ 0.7 eV, A = 6.5 
cV, B = 0.15 eV, C = 0.58 eV, F = 5 eV, F 2 = 6 eV) 
evidences that copper and oxygen Dzyaloshinsky vectors 



D o (0) = 



16 E t (jp x p y ) \e 



I 

= -2i{00\V so (O)\lz) =£([c.(dp) x c t (dp)] z + [c s (pd) x c t (pd)] z ), 

\ \/3 tpda . , . V3 tpda 
Cs,t{dPx) = 7T & s sm o ' CsAdPy) = 7T p \ cos o ' 

2 E s ,t{dp x ) 2 2 E St t{dp y ) 2 



D O (0) = 



3^2pi 



2p b pda 



Ldpx 



Idpy_\ 
e V ) 



3£,2pt 



?-V l pda 

~8 



sm 



cos" 



'-dpa i 
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can be of comparable magnitude. 



D. Dzyaloshinsky-Moriya coupling in La2CuC>4 

The DM coupling and magnetic anisotropy in 
La2CuC>4 and related compounds has attracted consid- 
erable attention in 90-ths (see, e.g., RefsQHH), and is 
still debated in the literatureii*^ In the low-temperature 
tetragonal (LTT) and orthorhombic (LTO) phases of 
La2CuC>4, the oxygen octahedra surrounding each copper 
ion rotate by a small tilting angle (Sltt ~ 3°, Slto ~ 5°) 
relative to their location in the high-temperature tetrag- 
onal phase. The structural distortion allows for the 
appearance of the antisymmetric Dzyaloshinsky-Moriya 
interaction. In terms of our choice for structural pa- 
rameters to describe the Cui-0-Cu2 bond we have for 
LTT phase: 6 — it; Si = 6 2 = § ± Sltt for bonds 
oriented perpendicular to the tilting plane, and 9 = 
±(tt — 2Sltt)',Si = 82 = f for bonds oriented parallel to 
the tilting plane. It means that all the local Dzyaloshin- 
sky vectors turn into zero for the former bonds, and turn 
out to be perpendicular to the tilting plane for the latter 
bonds. For LTO phased = ±(tt - V^<5lto);<5i = $2 = 
■| ± 8 lto- In such a case the largest (oc Slto) compo- 
nent of the local Dzyaloshinsky vectors (z-component in 
our notation) turn out to be oriented perpendicular to 
the Cui-0-Cu2 bond plane. Other two components of 
the local Dzyaloshinsky vectors are fairly small: that of 
perpendicular to Cu02 plane (y-component in our nota- 
tion) is of the order of S LTO , while that of oriented along 
CU1-CU2 bond axis (x-componcnts in our notation) is of 
the order of Sfr rn . 



IV. DM COUPLED CU1-O-CU2 BOND IN 
EXTERNAL FIELDS 

A. Uniform external magnetic field 

Application of an uniform external magnetic field hs 
will produce a staggered spin polarization in the antifer- 
romagnetically coupled CU1-CU2 pair 

<V 12 )=L = --^Dgxh s ] = X h s (25) 

with antisymmetric V r S'-susceptibility tensor: XaB = 

~ Xpa ■ One sees t na, t the sense of a staggered spin po- 
larization, or antiferromagnetic vector, depends on that 
of Dzyaloshinsky vector^ The VS coupling results in 
many interesting effects for the DM systems, in particu- 
lar, the "field- induced gap" phenomena in ID s=l/2 an- 
tiferromagnetic Heisenberg system with alternating DM 
coupling. 6 Approximately, the phenomenon is described 
by a so called staggered s=l/2 antiferromagnetic Heisen- 



berg model with the Hamiltonian 

H = J^(§, ■ s i+ i) - h u s lz - {-l) l h s s lx , (26) 

i 

which includes the effective uniform field h u and the in- 
duced staggered field h s oc h u perpendicular both to the 
applied uniform magnetic field and Dzyaloshinsky vector. 



B. Staggered external field 

Application of a staggered field h v for an antiferro- 
magnetically coupled CU1-CU2 pair will produce a local 
spin polarization both on copper and oxygen sites 

(S i ) = - r [D«xh^] = X (i)h v , (27) 

J 12 

which can be detected by different site-sensitive meth- 
ods including neutron diffraction and, mainly, by nu- 
clear magnetic resonance. It should be noted that SV- 
susceptibility tensor is the antisymmetric one: xfj = 
— XpX ■ Strictly speaking, both formulas H25|M27fl wor k 
well only in a paramagnetic regime and for relatively 
weak external fields. 

Above we addressed a single CU1-O-CU2 bond, where, 
despite a site location, the direction and magnitude 
of Dzyaloshinsky vector depends strongly on the bond 
strength and geometry. It is clear that a site rather 
than a bond location of DM vectors would result in a 
revisit of conventional symmetry considerations and of 
magnetic structure in weak ferro- and antiferromagnets. 
Interestingly that the expression (|27f) predicts the effects 
of a constructive or destructive (frustration) interference 
of copper spin polarizations in ID, 2D, and 3D lattices 
depending on the relative sign of Dzyaloshinsky vectors 
and staggered fields for nearest neighbours. It should 
be noted that with the destructive interference the lo- 
cal copper spin polarization may turn into zero and DM 
coupling will manifest itself only through the oxygen spin 
polarization. Another interesting manifestation of oxy- 
gen Dzyaloshinsky-Moriya antisymmetric exchange cou- 
pling concerns the edge-shared Cu02 chains (see Fig. 2), 
ubiquitous in many cuprates, where we deal with an ex- 
act compensation of copper contributions to Dzyaloshin- 
sky vectors and the possibility to observe the effects of 
uncompensated though oppositely directed local oxygen 
contributions. Applying the staggered field along chain 
direction (O^) we arrive in accordance with Exd. (|27|I at 
a weakly antiferromagnetic O y ordering of oxygen spin 
polarizations which magnitude is expected to be strongly 
enhanced due to usually small magnitudes of 90° sym- 
metric superexchange. It should be emphasized that the 
summary Dzyaloshinsky vector turns into zero hence in 
terms of a conventional approach to DM theory we miss 
the anomalous oxygen spin polarization effect. In this 
connection it is worth noting the neutron diffraction data 
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FIG. 2: The fragment of a typical edge-shared CV1O2 chain 
with copper and oxygen spin orientation under staggered field 
applied along x-direction. Note the antiparallel orientation of 
oxygen Dzyaloshinsky vectors. 



by Chung et alSL which unambiguously evidence the oxy- 
gen momentum formation and canting in edge shared 
Cu02 chain cuprate Li2Cu02- 



a conventional formula for a spin-dipole contribution to 
local field 



Hn — —9s[l>B ^ 



3(r l • Si)rj - rfsi 



and making use of an expression for appropriate matrix 
element 

(p l \ 3XaXf3 ^ r — \pj) = ~ (J^) (j>i\3lah-2Sa!3\Pj) ■ 



2p 



2/1 \ r 3 

2p 



~>\r 3 / ^2 SaiSf33 ' + 2 SajSpi ~ Saf}5i ^ ^ 28 ^ 



we present a local field on the 17 nucleus in CU1-O-CU2 
system as a sum of ferro- and antiferromagnetic contri- 
butions as follows 26 



H n = a" • (S) + A V • (V) 



(29) 



where 



S S S V .]/ y 

A =A (dp)+A (pd);A =A (pd) - A (dp), 



C. Intermediate oxygen 17 O Knight shift as an 
effective tool to inspect DM coupling for Cu-O-Cu 
bonds 

As an important by-product of the cuprate activity we 
arrived at a significant progress in different experimen- 
tal methods and new opportunities to elucidate subtle 
details of electron and spin structure. In particular, it 
concerns the oxygen 17 NMR-NQR as an unique local 
probe to study the charge and spin densities on oxygen 
sites^ In this connection we point to papers by R. Walst- 
edt et alm& as a first direct observation of anomalous oxy- 
gen hyperfine interactions in generic cuprate La2Cu04. 
With the approaching transition to the ordered magnetic 
phase, the authors observed anomalously large negative 
17 Knight shift for planar oxygens which anisotropy re- 
sembled that of weak ferromagnetism in this cuprate. 
The giant shift was observed only when external field 
was parallel to the local Cu-O-Cu bond axis (PL1 lines) 
or perpendicular to Cu02 plane. The effect was not 
observed for PL2 lines which correspond to oxygens in 
the local Cu-O-Cu bonds which axis is perpendicular 
to in-plane external field. The data were interpreted as 
an indication of oxygen spin polarization due to a local 
Dzyaloshinsky-Moriya antisymmetric exchange coupling. 
However, either interpretation of NMR-NQR data in such 
low-symmetry systems as La2Cu04 needs in a thorough 
analysis of transferred hyperfine interactions and a re- 
visit of some textbook results being typical for the model 
high-symmetry systems. First we draw attention to spin- 
dipole hyperfine interactions for O 2p-holes which are 
main participants of CU1-O-CU2 bonding. Starting from 



Af^dp) = A^[3c t (d Pl )c t (d Pj ) - \c t (d P )\ 2 S^] , 
Af^pd) = A^[3c t ( Pl d)c t ( Pj d) ~ |c t (pd)| 2 <y , 
AYjidp) = A^[3c s (dp7j^{dpj) - {c s (d P ) ■ c t {dp))5 l3 ] , 

AYjipd) = AWficifadfafaeC) - (c.(pd) • c t (pef))<%] , 

where A p ^ = |g s /is (pr) 2p ' ^ ne tilde points to a sym- 
metrization. Thus, along with a conventional textbook 
ferromagnetic (oc (S)) transferred hyperfine contribution 
to local field which simply mirrors a sum total of two 
Cu-0 bonds, we arrive at an additional unconventional 
antiferromagnetic difference (oc (V)) contribution which 
symmetry and magnitude strongly depend on the orien- 
tation of the oxygen crystal field axes and CU1-O-CU2 
bonding angle. With the CU1-O-CU2 geometry shown in 

Fig.l we arrive at a diagonal A tensor: 
A s xx = 2A p (3sin 2 |-1); A s yy = 2A p (3cos 2 A s zz 

and the only nonzero components of A tensor: 
AY.. = AY„ = 3A„sin( 

with 



2An 



xy yx 



A 



3 / tdpa_ \ , _ f ,0 

' _ 4 V e ' p _ p ' 
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where f„ is the parameter of a transferred spin density 
and we made use of a simple approximation E s t (dp x , y ) « 
e p . Generally speaking, we should take into account an 
additional contribution of magneto-dipole hyperfine in- 
teractions. 

The two-term structure of oxygen local field implies a 
two-term S-V structure of the 17 Knight shift 

^s^ss „v„i's 
17 K = A X +A X 

that points to Knight shift as an effective tool to in- 
spect both uniform and staggered spin polarization. The 
existence of antiferromagnetic term in oxygen hyperfine 
interactions yields a rather simple explanation of the 
17 Knight shift anomalies in La 2 Cu04 13 as a result 

of the external field induced staggered spin polarization 
„vs 

(V) = L = X Hext- Indeed, "our" local y axis for 
CU1-O-CU2 bond corresponds to the crystal tetragonal c- 
axis oriented perpendicular to Cu02 planes both in LTO 
and LTT phases of La2Cu04 while x-axis does to local 
Cu-O-Cu bond axis. It means that for the geometry of 
the experiment by Walstedt et alm& (the crystal is ori- 
ented so that the external uniform field is either || or 
_L to the local Cu-O-Cu bond axis) the antiferromag- 
netic contribution to 17 Knight shift will be observed 
only a) for oxygens in Cui-0-Cu2 bonds oriented along 
external field or b) for external field along tetragonal c- 
axis. Experimental data— agree with staggered magne- 
tization along the tetragonal c-axis in the former and 
along the rhombic c-axis (tetragonal a6-axis) in the lat- 
ter. Interestingly, the sizeable effect has been observed 
in La2Cu04 for temperatures T ~ 500 K that is essen- 
tially higher than T N » 300 K. Given L = 1, A p 0) » 100 
kG/spinji 3 . | sin 1 » 0.1, and f a m 20% we obtain « 6 
kG as a maximal value of a low-temperature antiferro- 
magnetic contribution to hyperfine field which is paral- 
lel to external magnetic field. This value agrees with a 
low-temperature extrapolation of high-temperature ex- 
perimental data by Walstedt et alm& Similar effect of 
anomalous 13 C Knight shift has recently been observed 
in copper pyrimidine dinitrate [CuPM(N03)2(H 2 0)2]n, a 
one-dimensional S=l/2 antiferromagnet with alternating 
local symmetry^ However, the authors did take into ac- 
count only the inter-site magneto-dipole contribution to 

A tensor that questions their quantitative conclusions 
regarding the " giant" spin canting. 

The ferro-antiferromagnetic S-V structure of local field 
on the nucleus of an intermediate oxygen ion in a Cui- 
O-CU2 triad points to 17 NMR as, probably, the only 
experimental technique to measure both the value, direc- 
tion, and the sense of Dzyaloshinsky vector. The latter 
possibility was realized earlier with 19 F NMR for weak 
ferromagnet FeF3iSi The experimental observation of an- 
tiferromagnetic contribution to the Knight shift provides 
probably the sole way to find out the problem of the sense 
of Dzyaloshinsky vector in cuprates. For instance, the 
negative sign of 17 Knight shift in L^CuO^i 3 - points 



^vs 

to a negative sign of X for Cui-0-Cu2 triad with 
A^y > 0, hence to a positive sense of z-component of the 
summary Dzyaloshinsky vector in Cui-0-Cu2 triad with 
geometry shown in Fig.l given 9 < it, Si — 62 ~ tt/2. 
This agrees with theoretical predictions both for cop- 
per and oxygen contributions based on Exps. l|17fl and 
(|'2ip. It should be emphasized that the above effect 
is determined by the summary Dzyaloshinsky vector in 
CU1-O-CU2 triad rather than by a local oxygen "weak- 
ferromagnetic" polarization as it was proposed by Wal- 
stedt et.alU 



V. SYMMETRIC SPIN ANISOTROPY 

Symmetric spin anisotropy for a CU1-CU2 pair is de- 
scribed by an effective Hamiltonian 

1 »^s\ 1 „ «V „ 
H an = ^SK S - -VK V (30) 

with kinematic relations Q. 

Depending on the sign of the anisotropy constants we 
arrive at two types of spin configurations minimizing the 
energy of spin anisotropy: the conventional twofold de- 
generated ferromagnetic state or unconventional multiple 
degenerated antiferromagnetic state. As a relevant illus- 
trative example one might refer to the axial anisotropy: 
H an = KS Z which stabilizes the |1 ± 1) doublet given 
K < or a set of superposition states 

* Q , = cosa|00) + e # sina|10) (31) 

given K > 0. The latter incorporates the limiting con- 
figurations 1 00) and 1 10), the Neel (-73GOO) ± |10)) and 

DM doublets (-73GOO) ± i|10)), and any other their su- 
perpositions. 

As usual, the term is processed in frames of a num- 
ber of simple model approximations. First, instead of 
the generalized form lj3U)l one addresses a pseudodipole 
anisotropy 

H an = S1K12 s 2 ■ (32) 

Second, one makes use of a trivial local MFA approach 
which applicability for s=l/2 spin systems is question- 
able. Indeed, calculating the classical energy of the pseu- 
dodipole two-ion anisotropy K(§i z )(s2 Z ) = jK((S z ) 2 — 
(V z ) 2 ) and respective quantum energy K{s\ z S2 Z ) = 
\K{{Sf) - (V 2 )) for a two-ion state cosa|00) + sinallrt) 
induced by a Neel-like staggered field || n we arrive at 
Edass = -\K sin 2a • n 2 z and E quant = \K{\ - 2 sin 2 a • 
n 2 .), respectively, that evidences the crucial importance of 
quantum effects when addressing the numerical aspect of 
spin anisotropy. Note that the mean value (V z ) reaches 
the maximum (=1) on a set of superposition states 1311 1 . 
while (Vz) 2 does on a single Neel state x I' Q=7r /4.0 = o- 
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A. Effective symmetric spin anisotropy due to DM 
interaction 

When one says about an effective spin anisotropy due 
to DM interaction one usually addresses a simple classical 
two-sublattice weak ferromagnet where the free energy 
has a minimum when both ferro- (oc (S)) and antiferro- 
magnetic (oc (V)) vectors, being perpendicular to each 
other, lie in the plane perpendicular to Dzyaloshinsky 
vector D. However, the issue is rather involved and ap- 
peared for a long time to be hotly debatedi2*ifl*22i211 In our 
opinion, first of all we should define what spin anisotropy 
is. Indeed, the description of any spin system implies the 
free energy $ depends on a set of vectorial order parame- 
ters (e.g.,(S),oc (V),oc (T)) under kinematic constraint, 
rather than a single magnetic moment as in a simple fer- 
romagnet, that can make the orientational dependence of 
<J> extremely involved. Such a situation needs in a careful 
analysis of respective spin Hamiltonian with a choice of 
proper approximations. 

Effective symmetric spin anisotropy due to DM inter- 
action can be easily derived as a second order perturba- 
tion correction due to DM coupling as follows: 

-^on M = PHdmRHdmP j 

where P is the projection operator projecting on the 
ground manifold and R = (1 — P)/(E — Ho). For an- 
tiferromagnetically coupled spin 1/2 pair H^ 1 ma y be 
written as follows: 



of an antiferromagnetically coupled spin 1/2 pair in a 
N eel-like staggered field h v || n: 



with AKVj = jjDiDj provided |D| < J. We see that in 
frames of a simple MFA approach this anisotropy stabi- 
lizes a Neel state with (V) _L D. However, in fact we deal 
with an MFA artefact. Indeed, let examine the second 
order perturbation correction to the ground state energy 



771 DM 

H>n„ — — 



ID-nP 



ID x n| 2 



4(£|, - E g ) A(E ± - Eg 



cos 2 a . 



(33) 



where Ej_ = J; Eu = J cos 2 a + h v sin 2a; E g = 
J sin 2 a — h v sin 2a. First term in (|33|) stabilizes n || D 
configuration while the second one does the n _L D 
configuration. Interestingly that (£j| — E g )cos 2 a = 
(E± — Eg), that is for any staggered field E^ 1 docs 
not depend on its orientation. In other words, at vari- 
ance with a simple MFA approach, the DM contribution 
to the energy of anisotropy for an exchange coupled spin 
1/2 pair in a staggered field turns into zero. The con- 
clusion proves to be correct in the limit of a zero field as 
well. 

B. Microscopies of symmetric spin anisotropy 

Anyway, the H^ 1 term has not to be included into 
an effective spin anisotropy Hamiltonian ©. As con- 
cerns the true symmetric two-ion spin anisotropy (pseu- 
dodipole, or exchange anisotropy), its magnitude can be 
obtained if we take into account all other effects quadratic 
in spin-orbital coupling. 

At variance with the effective DM spin-Hamiltonian 
the symmetric spin anisotropy evolves from the higher- 
order perturbation effects that makes its analysis even 
more involved and leads to many misleading estimations. 
Similarly to DM interaction we deal with two competing 
contributions. The first is derived as a lowest order con- 
tribution which does not take account of orbital fluctua- 
tions for Cui.2 3d states. To this end we take into account 
the effects of spin-orbital mixing for ground state singlet 
and triplet 101 configurations perturbed by covalent ef- 
fects. Assuming the validity of conventional perturbation 
series we arrive at a modified expression for respective 
functions as follows 



J 



*101;SA/ = $101;SM + E C {«}( 25+lr ) 

Mr 



{n}-TSM 



{{n};VS'M'\V ao \{n};TSM) 
S ^ T , E 2S , +lr ,({n}) - £ 2S+lro (101) 



{n}-X'S'M' 



(34) 



r 



and arrive at the expressions for the tensorial anisotropy parameters as follows: 
I 



K 



« = E C{n } ( 3 ri) W 3 r 2 ) 

{n}r 1 ,r a ,r' 



({n}; 3 r 1 |Af|{n}; 3 r)({n}; 3 r|Af|{n}; 3 r 2 ) 
E 3r ,({n})~Es ro (Wl) 



(35) 
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K v 



{n}r!,r 2 ,r' 



! ri)c {n} ( 3 r 2 ) 



Biiv({n})-^ ro (101) 



(36) 



We see that Kfj and KZ are determined by the triplet- 
triplet and singlet-triplet mixing, respectively. Interest- 
ingly that for nonzero orbital matrix elements in (|35[) and 
we find 

(M^rxlAfKn}; 3 ^) = (M^rxlAfKn}; 1 ^) 

hence Kf- = KY, if we will suppose that Ei T *({n}) = 
i?3 r /({ro}), that is neglect the singlet-triplet splitting for 
r' terms, or respective exchange effects. It should be 
noted that the contribution of two-hole on-site 200, 002, 
and 020 configurations in 1)35(1 and H36[l is actually re- 
lated to an on-site, or single ion spin anisotropy. Thus 
we conclude that, strictly speaking, a simple two-site 
pseudodipole form of symmetric anisotropy |QJ fails to 
capture correctly all the features of spin anisotropy in 
our three-center two-hole system. Firstly it concerns the 
quantitative predictions and estimations. The contri- 
bution of two-hole two-site 110 and 011 configurations 
to spin anisotropy turns out to be nonzero, only if pd- 
exchange is taken into account. 

The second contribution is associated with orbital fluc- 
tuations for Cui, 2 3d states within a ground state 101 
configuration, and evolves from a third order combined 
effect of Cux,2 spin-orbital V so (Cui^) and effective Cui- 
Cu ? ex change coupling (see, e.g., a detailed analysis in 
Ref l28T) . Thus, we see that any decisive conclusions re- 
garding the quantitative estimations of symmetric spin 
anisotropy imply a thorough analysis of numerous com- 
peting contributions. 

VI. CONCLUSIONS 

In conclusion, we have revisited a problem of 
Dzyaloshinsky-Moriya antisymmetric exchange coupling 
in cuprates specifying the local spin-orbital contributions 
to Dzyaloshinsky vector focusing on the oxygen term. We 



have applied a scheme that provides an optimal way to 
account for intra-atomic electron correlations and sepa- 
rate the local contributions to Dzyaloshinsky vector. The 
Dzyaloshinsky vector and respective weak ferromagnetic 
momentum is shown to be a superposition of compara- 
ble and, sometimes, competing local Cu and O contri- 
butions. The intermediate oxygen 17 Knight shift is 
shown to be an effective tool to inspect the effects of 
Dzyaloshinsky-Moriya coupling in an external magnetic 
field. Anisotropic antiferromagnetic contribution to 17 K 
explains the anomalies observed in I^CuC^^ We pre- 
dict the effect of strong oxygen weak antiferromagnetism 
in edge-shared CuC>2 chains due to uncompensated oxy- 
gen Dzyaloshinsky vectors. Its experimental observation 
could provide a direct evidence of the oxygen DM cou- 
pling. We revisit the effects of symmetric spin anisotropy, 
in particular, those directly induced by Dzyaloshinsky- 
Moriya coupling. 

Finally it should be noted that the anionic contribution 
to Dzyaloshinsky vector is crucial for the very existence 
of the DM coupling in the pair of rare-earth ions (see, 
e.g., Yb 3+ -As 4 "-Yb 3+ triads in Yb 4 As3^i) because very 
strong spin-orbital coupling for rare-earth ions is diago- 
nalized within a ground state multiplet. 
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